PATH COUNT ASYMPTOTICS AND STIRLING NUMBERS 
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Abstract. We obtain formulas for the growth rate of the numbers of certain 
paths in a multi-dimensional analogue of the Eulcrian graph. Corollaries are 
new identities relating Stirling numbers of the first and second kinds. 



1. Introduction 

A graded, infinite, directed graph defines a dynamical system called a Bratteli- 
Vcrshik or adic system (see |Sj and the references given there for background on 
adic systems). One of the main problems concerning these systems is to determine 
their ergodic invariant measures, and the basic tool for attacking this problem is 
the asymptotics of numbers of paths in the graph with certain properties. 

A very interesting Bratteli-Vershik dynamical system is associated to the Euler- 
ian two-dimensional graph with vertices at the points of Z> . The paths in the 
Eulerian graph may be counted in many different combinatorial ways, but all at- 
tempts to extend those counts to the multi-dimensional analogues of the Eulerian 
graph have failed up to this moment; techniques applied to the two-dimensional 
Eulerian graph, such as the recurrence of random walks in two dimensions in pQ, 
the coding of paths by permutations in [5] , or the left-to-right monotonicity in [7] , 
are not available in dimensions higher than two. 

In this paper we make a first step in the multi-dimensional direction. We consider 
a multi-dimensional analogue of the Eulcrian graph and determine the asymptotics 
of the numbers of paths with certain properties. We plan to apply these asymptotics 
to studying the ergodic measures on these and related systems in future work. As 
a rather unexpected corollary of our asymptotics we get a family of new identities 
between Stirling numbers of the first and second kinds, which are classical objects, 
arise in a variety of combinatorics problems, and for which many general relations 
have long been known (see for example [91IT0]). 
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Now we formulate the results of this paper. Denote by e\ = (1, 0, . . . , 0), . . . , 
e„+i = (0, ...,0, 1) e the ordinary basis vectors in E n+1 . Let c e Z™^ 1 . 

We define the graph E c to be the graph with vertices i S Z™^ 1 and the following 
directed multiple edges. Each vertex i is connected just ton + 1 vertices: to vertex 

i + ej by Cj + i n+1 edges if j < n + 1 and to i + e n+ i by c+i + ii + «2 H h *n 

edges. 

Such a graph models a random walk on the graph G n +i consisting of n + 1 loops 
based at a single vertex, with a special kind of opposite reinforcement: when the 
walker takes a step in the direction of e n +i , the numbers of edges in all the other 
directions (ej for j = 1, . . . , n) are incremented by 1; but when the walker takes a 
step in the direction of ej for some j = 1, . . . , n, only the number of edges in the 
direction of e n+ i is incremented by 1 (cf. [6]). When n = 1 and c = (1, 1), E c is 
the standard two-dimensional Eulerian graph. The corresponding two-dimensional 
Eulerian adic system was introduced in [T| and studied further in j5] and [8]; see 
also [HE]. 

For each i G ^"g 1 denote by A c (i) the number of paths in the graph E c from 
to i. The standard Eulerian numbers (see [3j|8]) are the numbers A(i\,i2) = 
-^(0,0) (*i 1*2) with n = 1. In this paper, we give a formula for the leading term of 
the growth rate of A c (i\, . . . , i n , i n +i) as i n +i tends to infinity. We show that if 
m = i\ + • • • + i n , then 

A c (i u . . .,i n ,i n+1 ) ~ B c (ii, . . .,i n )(c n +i + m) ln+1 

as i n +i — > 00 and 

(1-1) £ Cl ,... )Cn+I (ii,...,i n ) = " +1 m , 51 II ( c /(i) + c »+i +i _ !) » 

f£F j=l 

where F denotes the set of all maps / : {1, . . . , m] — > {1, . . . , n} such that |/ _1 (j)| = 
ij for j = l,...,n. 

As a corollary we obtain new identities relating Stirling numbers of the first and 
second kinds; for 1 < k < n, < r < k, we get 



r + n — k — 1 

r 



s\(n, r + n — k) 



(1.2) J^fm + n-k\^fi + Ti-k + m-l\(-iy n+r -\ 

1 JM i ) Ti^+i — * + 1)1 X 



?n— v 7 i— 

S2(m + 1, r — i + l)si(n, i + n — k + m), 

where s\,S2 are Stirling numbers of the first and second kinds: see Theorem 13.101 
and Corollary 13. Ill 

Acknowledgment. We thank R. Stanley for telling us about theorem of Frobenius 
(Theorem 13. 6|) and the referee for suggestions leading to the improvement of the 
exposition. 
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2. Path counts and limits 



2.1. Finding the leading term. The numbers A c (i) of paths in the graph E c 
from to i satisfy the recurrence relation 

n 

(2.1) A c (i) = (c„+i + i\ H h i„)A c (i - e„+i) + E( c 3 + i n +i)A:(i - e^). 

3=1 



We are interested in the limit 

A c (ii, . . . ,i n ,i n +i) 



(2.2) S^,...,^) = _ lim 



i„+i-^oo (c„ + i + n) l "+ 1 
Proposition 2.1. Lei it = (1, . . . , 1) G ZtJ. . TTiera 



(2.3) B c (u) = E E II ( c e(^) + 3i + • ■ ■ + 3k) 

ees„3i,-,3n=o fe=i 



c ra+i + k - 1 

Cn+l + n 



Proof. To form a path from to i = (1, . . . , 1, i n +i) in E c , we need to take one step 
in the direction of each e^, k = 1, . . . ,n, and i n +\ steps in the direction of e n +\- 
Thus each such path will correspond to a permutation £ of 1, 2, . . . , n which will 
determine when the step in the direction of each e&, k = 1, . . . ,n, is taken, together 
with a partition i n +\ = Ji + ■ • • + J n +Jn+i of i„+i into nonnegative integers, which 
will determine how many consecutive steps are taken in the e n +i direction between 
steps in other directions. Each time that j steps are taken in the e„+i direction, 
all the cfc, 1 < k < n, are incremented by 1; and each time that a step is taken in a 
direction other than e n +i, c n +i is incremented by 1. Therefore, 
(2.4) 

A c (l,...,l,in+i) 

= E E C n+l( C «(l) + Jl)( C "+l + 1 ) J2 ( C «(2) +J1 +h) ■ ■ ■ 

£6S„ 3iH h3„ + i=i„ + i 

• • ' (Cg(„) +.71 H h Jn)(c„+1 + + 1 



E E 

£eS„ 3iH — hi„+i=i„+i 



(c„+i 
(c„+i -+ 



J] (c„+i + k - l) jk (c m +;)! + ■■■+ jk) 
(Cn+l + k-l) 



.fe=l 



(c„+i + n) J " +1 



)-e e n 

in+i n 

r- +1 E E n 



js>0 k = l 

Jl + ---+Jn=0 



Cn+l + n 

(c n+1 + k-iy Jk 



Cn+i + n 



(Cf(jfe) + jiH h jfe) 



(c£(k) +ii + hife). 



This reasoning gives Formula (|2.3|) . Clearly, B c {u) is a convergent sum. □ 
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Similarly, existence of the limit in Formula (|2.2j) for any (ii, . . . , i n ) can be shown 
as follows. Let m, F, be the same as in Section [T] Then 

Ac[j>li ■•■i'iii in+l ) = 

= H H 4+i( c /(i) + ji)( c «+i + !) J2 ( c /(2) +ii +J2) • • • 

f&F jH hj m +i=i„+i 

• ' ' (c/( TO ) + ii H h jm)(c„+i + m) Jm+1 

= (c„ +1 +m)^£ ]T n 

f£F 3a>0 fe=l 

>'i+-+i»=o 

Dividing by (c n +i + m) 4n+1 and letting i„+i —> 00 again leaves a sum of convergent 
series. Thus we get 

00 n 

(2.6) B c (ii,...,i n ) = 5Z II ( c /(fe)+ii + --'+ifc) 

feF h,...,j n =o k=i 



(c 



n+l 



1) 



c„+i + m 



( c /( 



fc) + h + 



+ Jk)- 



/ Cn+1 + k-i y 

\ c n+1 +n J 



Wc shall give two formulas for B c {u). In this section we evaluate B c (u) using 
a two-dimensional formula from [8]. In Section 13.11 another formula is obtained 
by using derivatives of geometric series. Comparing the two formulas we get new 
identities involving Stirling numbers of the first and second kinds (Corollary 13. 11 1) . 

Theorem 2.2. We have 

(2.7) B c {u) = (Cn+1 | nJ £ 1] ( c «^) + c "+i + k !) • 

' £es„ fc=i 



2.2. Proof of Theorem 12.21 Denote by erj(c), i = 1, . . . , n, the elementary sym- 
metric functions in c, 

(2.8) ai(c)= X ^••• c i i - 

l<3i< — <J»<n 

Denote cro( c ) = 1- Clearly, the function B c (u) is a symmetric polynomial in c of 
the form 

n 

B c (u) = ^ i] -( n ~ *) ! °~i{c)<Xi,n{cn+i) = n\ <r n (c)a n ,„ (c„+i ) + 

(n - 1)! <7 n _i(c)a„_i. n (c„ + i) + ... + n!cr (c) ^0,71 (^n+l) 1 
where a^ n (c n +i), i = 0, . . . ,n, are suitable functions in c n +i. For example, 

c„+i + fc-l v 



a ,n(c„+i) = ^ JJ O'i H hjfc 

(2.10) Jl '- J " =0 fc=1 

/ \ TT I c n+l ' 

a n ,n(c„+i) = 2^ 

ii,-,j„=o fc=i 

and so on. 



Cn+1 + n 



Cn+1 + n 
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Example 2.3. For n = 3, 

OO 

B (c u c2,c 3 ){u) = ^2 I( c l + Jl)( c 2 +J1 + J2) + (c 2 + jl)(ci +.n +j 2 )] X 

<>•"> x7^y fail)* 

= 2! cic 2 o 2 , 2 (c 3 ) + (ci + c 2 ) ai, 2 (c 3 ) + 2!a , 2 (c 3 ), 



where 



« 2 , 2 (c 3 ) = 

3ijja=0 



C3 \ J VC3 + 1 XJ2 



(2.12) ai, 2 (c 3 ) = ( 2 h+h) 



c 3 V 1 (c 3 + l^ 32 



c a + 2J \c 3 + 2 



31,32=0 

«0, 2 (C3)=^E o A(i 1 +i 2 )(^) (J^) 



Our goal is to calculate all the coefficients oci, n (cn+i) an d show that the result 
gives Formula (|2.7[) . To do so let us calculate the number 



jk 



(2.13) B CuC2 (n) = II {^+h + ---+3k)[^7 

h,...,j n =o fc=i v c 2 +n y 

The function 5 Cl , C2 (n) has the form: 

(2.14) B Cl , C2 (n) = c"a ni „(c 2 ) + c"^ 1 a„_i ! „(c 2 ) H ha ,n(c 2 ), 

where the functions aj,„, i = 0, . . . ,n, are the same functions of one variable as in 
()2.9|) . (To see this, note the identical form of the two products in formulas (|2.3j) 
and (|2.13j) . and note that the sum over the symmetric group is absorbed into the 
coefficients that do not involve c n +i-) 



C2 yvc 2 + i xJ2 



Example 2.4. We have 

oo 

B CUC2 (2)= £ (ci+h)(ci+ji+32 / , 
(2.15) h ±L Vc 2 +2y \c 2 + 2 

= c\a 2a {c 2 ) + ci 0:1,2(02) + 0-0,2(02), 
where the functions ao, 2 , ai, 2 , a 2 , 2 are the same functions of one variable as in 

mm). 



To calculate the functions Oi jrl and prove Theorem 12.21 it is enough to calculate 
B cuc 2 {n). 



Lemma 2.5. We have 

fe=i 



(2.16) B CuC2 (n) = {C2 + , n TT ( Cl + c 2 + fc-l) 
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Proof. Let us use formula (1.1) of [8] to evaluate 

(2.17) U^.taT^- 

i-s-oo [c-2 + n) 1 

Formula (1.1) for A Cl ^ C2 (n, i) has n + 1 summands and the leading summand is 

(2.18) (* + *+ n_1 ) («* + »)** 

This term after division by (c 2 + n) t gives the right hand side of (|2.17|) . The lemma 
is proved. □ 



Lemma [2751 gives Theorem 12.21 

2.3. Formula for B c (i). Formula (|2.7[) for S Cli ... )Cn+1 (l, . . . , 1) easily gives a for- 
mula for B Cl ,..., c „ +1 (ii, • ■ • , V) for an y •••)*'«) • 

Proposition 2.6. Let m ; F, be the same as in Section [3 T/iera 

(2.i9) ^....^(ix,...,^) = {cn+ \ir e n (c/o)+c n+1+J -i) . 



3. Stirling numbers and related polynomials 

Theorem 12.21 has as corollaries interesting identities involving Stirling numbers 
of the first and second kinds. 



3.1. Theorem 12.21 in new notation. Denote 

q — a d 



-"LI 5 LI — 'LL r\ 

q du a 

Then Theorem 12.21 says that for arbitrary numbers c\ , . . . , c„ we have 
(3.1) 

(n \ n _^ n n 

£ T[^ia) + Da + D a+1 + --- + D n ) m-__ = n( c -H+9-«)- 

creS,, a=l / 6=1 ° ' crGS„ q=1 



Example 3.1. For n = 2, Formula (|3.1j) says 

[(ci + £ 2 )(c 2 + D 2 + Di) + (c 2 + D a )(ci + D 2 + £>i)] : 



1 — Ml 1 — U 2 

a 2 

Y [( Cl + g - l)(c 2 + 9 - 2) + (c 2 + 9 - l)(ci + 9 - 2)] 
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1 




1 


1 


- 1*1 


1 


- U 2 




1 




1 


1 


- 1*1 


1 


- u 2 




1 




1 


1 


- Ul 


1 


- 1*2 



This statement is equivalent to the identities 

t 

2 ' 

= q - (q-2 + q-l) , 

D 2 (D 2 + D 1 ) - : = 2L(g-l)(g- 2 ). 

1 — 1*1 1 — 1*2 2 

3.2. Formula for (i*<9 u ) fe (l — it) -1 . We make some observations that are useful for 
handling expressions such as the left side of Formula (|3.ip . 

Let q be a variable, a a positive integer. Set 

(3.2) it = , d u = — . 

q ou 

Let A(i, j) = Ao.o(i,j) be the standard Eulerian numbers as in Section [T] 
Theorem 3.2. For any positive integer k we have 

(3.3) Wil-u)- 1 = a-"- 1 J2 A(i,j)q i+1 (q-ay + 1 . 

i+j—k-i 

Example 3.3. 

ud u {l - uy 1 = a~ 2 q(q~a), 

(3.4) (iicy 2 (l - u)- 1 = a' 3 {q 2 (q - a) + q(q - a) 2 ) , 

(ucy 3 (l - u)- 1 = a" 4 (q 3 (q - a) + 4q 2 (q - a) 2 + q(q - a) 3 ) . 



Theorem 13.21 follows by induction from the following lemma. 
Lemma 3.4. We have 

d q(q — a) d 

u- 



du a dq 

Proof of Theorem \3.2l 



* i+j=k-l 

- 2 J2 A(i,j)(i + l)q i+1 (q-ay+ 2 



a 

i+j=k-l 

- fc - 2 A(i ) i)(j + l)q i+2 (q-ay +1 



a 

fe-2 Af„ „ i \ /„ _, i \ „r+l/„ „\s+l 



a 

r-\-s=k 

,-fe-2 



£ A(r,«-l)(r + l)(f +1 («-a)' 

+s=fc 

£ i4(r-l,*)(* + l)g r+1 (9-a)* +1 
£ A(r,s) q r+1 (q-aY 



a 

r-\-s=k 

-fe-2 m„ „\„r+ir„ ^h-i 



r+s— fc 
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T k (q,n) = ^(-l) fe * S 2(M* 



□ 



3.3. T-polynomials. We shall reformulate Theorem 13.21 in terms of a slight vari- 
ation of Eulerian polynomials (see [3j p. 244]), which we call T-polynomials. 

Let s 2 (fc, m) denote Stirling numbers of the second kind, which count the number 
of partitions of a fc-element set into m nonempty subsets, so that 

sa(M) = l, 
s 2 (2,l) = l, s 2 (2,2) = l, 
s 2 (3,l) = l, s 2 (3,2) = 3, s 2 (3,3) = l. 

We have 

s 2 (fc, m) = s 2 (fc — 1, m — 1) + ms 2 (fc — 1, m) . 
We will be interested in the polynomials 

k 
i=l 

Example 3.5. 

T x {q,n) = l, T 2 (q, n) = 2q - n, T 3 (q, n) = 6q 2 - 6qn + n 2 . 

Theorem 3.6 (Frobenius, [3J p. 244]). We have 

T k (q,n) = Mhj)q l (q-ny . 

i+j=fc-i 

3.4. A-polynomials. We will also need A-polynomials. 

Looking at the right side of Formula p.l[) suggests investigating related polyno- 
mials and Stirling numbers of the first kind. 

Stirling numbers of the first kind are the coefficients in the expansion 

n 

(?)n = ^2s 1 (n,k)q k , 

k=Q 

where (q) n is the falling factorial 

(q) n = q(q-l)(q-2)---(q-n + l). 

We have 

si(n,fc)= {-l) n ~ k i 1 i 2 ---i n -k ■ 

l<ii<---<i n _k<.n— 1 

For natural numbers 1 < k < n, define a polynomial 

A„, fc (q) = ^ (q-ii)(q-i2)---(q-i n -k)- 

l<ii<---<? ra _fc<n — 1 

Define A 7hn (q) = 1. We have A nt x(q) = (q — 1) . . . (q — n + 1) and 
An,fc(g) = E (*"^ +m )*l(n,* + ro)<r. 



m— 
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Example 3.7. 

A 2 , 2 (g) = l, A 2)1 (g) = g-1, 
A 3 , 3 (g) = l, A 3 , 2 (g) = 2q - 3, A 3 ,i(g) = g 2 - 3g + 2. 

3.5. Relation between V and A polynomials. 

Example 3.8. We have the following identity: 

( i- 1 ®?^) o \ / A22(?) x _ / A22{q) 

y a-»g)r a («, 2 ) a- 1 ©^^) j ^ A 2|1 (g) ^ ^ A v (g) 

In other words, we have 

'2\ r 1 (g,2)A 2i2 (q) 



1J 2 
2\r 2 (g,2)A 2 , 2 (g) , A\r 1 (g,2)A a ,i(g) 



A 2 , 2 (q) , 

2J \ij =A 2 ,i(g). 

Example 3.9. We have the following identity 

/ 3_1 (?) rife3) \ / A 3>3 (g) \ / A 3 , 3 (g) 

3- 2 ©r 2(9 ,3) 3-^(1)^,3) o A 3)2 (g) = A 3 , 2 (g) 

^ 3- 3 (^)r 3 ( 9 ,3) 3- 2 (^)r 2 ( 9 , 3 ) 3-1(1)^,3) J \ A 3)1 (g) / \ A 3il (g) 

In other words, we have 

3^r 1 (g,3)A 3 , 3 (g) 

~ ^3,3 (.gj , 

A 3 , 2 (g) , 



IJ 3 
3\ r 2 (g,3)A 3i3 (g) , /2\ T 1 {q, 3) A 3)2 (g) 



3 2 Vl 



3^ r 3 (g,3)A 3 , 3 (g) ^ r 2 (g,3)A 3 , 2 (g) /l\ I\(g, 3) A 3 ,i(g) . 

3 J 3* + v 2 ) 3^— + [l ) 3 = AsAq) 



More generally, we have the following relation. 
Theorem 3.10. For any 1 < k < n we have 
(3.5) A n , fe (g) = ]T ^ + ™W^n)A„, fe+m (g) 



m=0 



,1+m 



Proof. Recall the notations and defined at the beginning of Section 13.11 

Denote 

(3.6) 

nn = n i n( c J)= e i](a j +a j+ i+--+a 1 )u„. 

fc=l b l<h< - <i k <n \j=l J 

For example, 

/i(g) = (nD n + (n - l)-D n -i + • • • + Di)n n . 
Denote ct;(c), i = 1, . . . , n, the elementary symmetric functions in c = (ci, . . . , c„) 
as before. 
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The left hand side of Formula (|3.1|) equals 

n 

]T k\(n-k)\a n - k (c)I*(q). 



fc=0 

The right hand side equals 



— V] k\(n - k)\ (T n _ fe (c)A„ + i n _ fe+ i(g) 



By comparing coefficients of the symmetric polynomials in these two equations we 
get 

q n 

I n (o) = r A„+i !n _ fe+ i(g) . 
n! 

We have 

A„ + i j „_ fc+ i(g) - A nj „_ fe (g) = (<? - n)A„ :)l _ fc +i(g) 

and 

*.>-§(::})«-.<«> K^)- 

Hence 



fc-i 



( 3. 7 ) /jw-^-w^-eC:^.*" H^t^ 



The left hand side of JO is 



3=0 



I'M - ^-i(9)rr- = - 'n-iteP 



— A n+ i n _fc +i (g) - — - — — A n n -k(q) — 

77 i in — 111 71 



(n-1) 

,n-fe+l(s) - 



— (A„+i. n _ fc+ i(g) - A„.„_ fc (<7)) 
n! 



The right hand side of (|3.7[) equals 

fc-i 



(9). 



fe-i 



§(::^-<"K't^) 



3=0 



These two equations prove Theorem 13. 101 □ 



Since A„ jn (g) = 1 and Ti(q, n) = 1 for all n, identity (|3.5|) allows us to express 
the polynomials A n ^k(q) in terms of polynomials Tk(q, n) and vice versa — see the 
examples above. 
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3.6. Identities involving Stirling numbers. The conclusion of Theorem 13.101 
can be written in the following form: 



1 > /to + k\ 1 

(3.8) A n ,k(q) = r > . 1 — r m+ i 5,n A n ,i; +m g , fc = l,...,n. 

n — k \ m + 1 / n m 

m— 1 x 7 

Comparing coefficients of g r of the two sides yields identities involving the Stirling 
numbers si(n, k) and S2(n, fc) of the first and second kinds, indexed with 1 < k < n 
(and equalling outside this range). 

Corollary 3.11. For 1 < k < n, < r < k, we have 

W + n — k — I s 



r 



Si(n, r + n — k) 



(3.9) A (m + n - k\ ^ (i + n -k + m-l\ (-1)"**-* 

m=0 v 7 i=0 v 7 

S2(m + 1, r — i + l)s±(n, i + n — k + m). 



For r = 0, this says 

,. . A (m + n-k\ (-l) m . , . 

n — K) — si(n,n — K)+ > sun, n — k + to), 

n V to + 1 / n 

m— 1 v ' 

and hence 

(3.10) Sl (n,n-fc) = - V ( m + n ~ k \-l) m Sl (n,n-k + m), k = l,...,n. 
k \ m + 1 I 

m—l v ' 

This contains the known formulas 

s\{n,n — 1) = — ( J and si(n, n — 2) = — n{n — l)(n — 2)(3n — 1). 
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